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1. INTRODUCTION 


The purpose of this appendix is to explain the method used in computing the 
structural influence coefficient matrix (or [C] matrix) of the computer program of 
Reference 1 (Appendix A of the Summary Report). This [Cl matrix is computed for 
complete wing-body-tail configurations as shown in Figure 1 by assuming that all major 
airplane components can be structurally represented by a slender beam called the elastic 
axis. Referring to Figure 1, a structural influence coefficient (i.e. any element Cjk of 
matrix [Cl ) is defined as the rotation about the Y-stability axis at panel "j" induced by 
a unit load on panel "k". Symbols used in this appendix are defined in Section 2. A 
description of how a structural breakdown is performed in detail is presented in Section 3. 
The method used in computing the elements of the structural influence coefficient matrix, 
[C] , is developed in Section 4. Validation of this method is discussed in Section 5. 





(b) Definition of C]|< 

Figure 1 Definition of Elements of the Structural Influence Coefficient Matrix 


2. SYMBOLS 


The 

units used for the physical quantities defined in this 

paper are given both in the 

Internationa 

1 System of Units (SI) and the U.S. Customary Units, 

• 

Symbol 

Definition 

I Dimension 

c ] 

structural influence coefficient matrix 

rad ./lb 
I (rad ./N . ) 


structural influence coefficient (change in 
of attack of panel j due to a unit load on 

angle 


panel k ) 

rad ./lb. 

— 


(rad ./N .) 

c 

s 

c.p. 

structural root chord 
panel center of pressure 

in.(m) 

El 

bending stiffness 

lb. -in. ^ 



(N.-m 2 ) 

GJ' 

torsional stiffness 

Ib.-in. 2 
(N ,-m 2 ) 

L i 

length of i-th elastic axis segment 

in . (m) 

M 

moment 

in-lb.(m.-N.) 

P 

force 

lb. (N.) 

T 

torque 

in-lb.(m.-N.) 


“ total streamwise rotation (angle of attack) 

of a point on the elastic axis rad. 

change in streamwise angle of attack 

of an elastic axis segment rad. 



Symbol 

Definition 

Dimension 


sweep angle of i-fh elastic axis segment 

rad. 

9 n 

bending angle due to a bending moment 

rad. 

*t 

bending angle due to a torque 

rad. 


twist angle due to a bending moment 

rad. 

0 T 

twist angle due to a torque 

rad. 



3. DESCRIPTION OF STRUCTURAL BREAKDOWN 


For purposes of analysis, the entire internal structure of the airplane is reduced to 
a number of elastic axes. For example, the elastic axis of the wing becomes a cantilever 
beam and represents the spars and other internal structural components. At each incremen- 
tal distance outboard of the aircraft centerline, the elastic axis assumes the equivalent El 
and GJ values for the particular components that the distance represents (see Figure 2a). 

The elastic axis is broken up into small segments called link lengths and in the 
analysis each link length is considered to be a cantilever beam (see Figure 2b). Each 
panel on the wing is assumed to be assigned to some segment of the elastic axis in a 
direction perpendicular to the elastic axis (see Figures 2b and 3). Attachment of the panel 
to a link length is accomplished by assuming an infinitely rigid arm connecting the panel 
center of pressure and the endpoint of the elastic segment. The term infinitely rigid, 
as used here, means that whatever rotations are induced in the elastic axis segment to 
which the jth panel is assigned, the jth panel also undergoes these same rotations. This 
infinitely rigid arm serves as a means to translate the load to the outboard point i on its 
attached elastic axis segment (see Figure 2b). 

Calculation of the structural influence coefficients is facilitated by the use of 
influence coefficients internal to the program. These internal influence coefficients are 
called Elastic Axis Flexibility Vectors. They are found as follows. First, a unit force, 
unit moment and unit torque (with positive orientation) are placed at the outboard end of 
the ith link length (see Figure4). Second, with this unit loading condition, the rotations 
induced at the loaded point i are calculated by considering the variance in bending and 
torsional stiffness inboard of the loaded point. This procedure is repeated for every link 
length starting at the outboard end of the first segment and proceeding out the axis until 
the end is reached. This procedure is discussed in more detail in Section 4.. 

Every effort was made to make the program as general as possible. As a result, the 
position of an aerodynamic panel can be anywhere within a 360° arc from its assumed 
point of attachment on the elastic axis (see Figure 5a). 

The center of pressure of the panel may have the same coordinates as its attachment 
point on the elastic axis if this is necessary (see Figure 5b). Discontinuities in the elastic 
axis can also be tolerated providing they are handled as the end of an elastic axis link 
length (see Figure 6). The link lengths may be as long or as short as desired and these 
lengths may vary from link to link. 

A positive loading condition is illustrated in Figure 4. A positive force is up as a 
lift force. Torque is positive when it induces a positive angle of attack for the wing, and 
a positive moment tends to bend the wing upward. 
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a) Shaded Area Shows Region Where Equivalent El and GJ Are Found. 



b) Panel Connection to Elastic Axis Link Length and Positive Loading Condition 


Figure 2 

. Example of Elastic Axis, Link Length and Panel Arrangement. 




Figure 3 Examples of Panel Assignment to the Elastic Axis 



Figure 4 Description of Unit Loading Condition 



(a) Aerodynamic Panel j May Be Oriented Anywhere on the 
Arbitrary Arc Shown. 


» Y 



(b) Coordinates of the Panel Control Point May Be the Same as the 
Coordinate of the End of an Axis Segment. 

Figure 5 Panel Assignment Possibilities 



Figure 6 Discontinuities to be Handled as Segment End Points 





4. METHOD OF ANALYSIS 


This section gives a discussion of the step-by-step procedures used in calculating 
the structural influence coefficient matrix of wing-body-tail configuration. The procedure 
for the wing alone case is described in Paragraph 4.1 . The procedure for the complete 
airplane case is described in Paragraph 4.2. 

4.1 Method for Computing the Stru ctural Influence Coefficient Matrix for th e Wing 

Calculation of the length of each elastic axis segment is the first computation com- 
pleted by the program. This is done by considering each segment of the axis to be the 
hypotenuse of a right triangle. Then using the coordinates of each end of the segment 
whose length is to be calculated, the following relation is used: 

j L i-1 = [*V X i-i> 2 + (y? - Xr-i) 2 p 0) 

I 

The coordinates used in Equation (1) are defined in Figure 7. 



Figure 7 Unit Loading System on Elastic Axis 

As stated in Section 3, calculation of the elastic axis flexibility vectors facilitates 
the computation of the structural influence coefficients for the aerodynamic panels. 

The procedure involves the placing of a unit force, unit moment and unit torque (with posi- 
tive orientation) at the outboard end of the i-th segment of the elastic axis (see Figure 7). 
The sweep angle for each elastic axis segment is found using the formula: 

■ y ;-i =,an "' ] • (2) 
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The unit torque-T applied at i is transferred in much the same way as the moment, 
giving the following equations for the bending angle dj and the twist angle, 0J : 



(sin ;_]) 0-) 


El. 

J 

(cosy .,) (L) 

Sr * 1 


GJ. 

J 


(7a) 


(7b) 


The next step in the computer program is to project the angles of twist and bending 
into the streamwise direction, giving the change in streamwise rotation, A a at the end- 


point j + 1 of elastic axis segment j. 



A “iyf'V c “ Y rV !n Y j 

(for unit force P) 

(8a) 

‘“My." -M cos Y j'«M sin Y j 

(for unit moment M) 

(8b) 

^ a Tj + | = 4>t cos Yj — 0 t s * n Y j 

(for unit torque T) 

(8c) 


The total streamwise rotation at any point on the elastic ax.i$ is simply the summation 
of the changes in rotation of all the points inboard of that point added to the change in 
rotation of that point. For example: 




Vm"^ 1 ; 

-(for unit force P at point i) 

(9a) 




J+ * / ' n 

(For unit moment M at point i) 

(9b) 

a = j+1 

Vi r2l da T n 

(for unit torque T at point i) 

(9c) 



Equations 9a, 9b, and 9c represent the desired elastic axis flexibility vectors. 
These equations hold as long as j + 1 is less than or equal to i; that is, for all segments 
inboard of the loading point i. If j + i is greater than i, then point j + 1 is outboard of 
point i, and the streamwise rotation of a point outboard of the loading point i is the same 
as that of point i . 

Now that the elastic axis flexibility vectors (Equations 9) have been calculated, 
the structural influence coefficients can be calculated. First, an IASIGN vector is con- 
structed which assigns each aerodynamic panel to an endpoint of the elastic axis (see 
Figure 11). 



Figure 11 Panel Attachment to Assigned Point on the Elastic Axis 


Construction of IASIGN involves judgment of how the actual structure should be interpreted 
in terms of the disjointed elastic axis and aerodynamic panel model which is used here. A 
unit load applied on any such panel K can be transferred to its assigned elastic axis end- 
point i by finding the torque arm T and moment arm M as follows: 


T = (y^ - /j) sin 7 ; _ 1 - (x k - x.) cos 7 M 


(10a) 


M = (y|< - y;) cos + ( x k" x j) sin 7 j_j 


(1 Ob) 


The streamwise rotation of the elastic axis can now be found by a linear combination 
of the elastic axis flexibility vectors of Equations (9). The rotation of any panel j due to 
the unit load on panel k is the same as the rotation of the endpoint of the elastic axis to 
which panel j is assigned. Therefore, if panel j is assigned to point j + 1 and panel k to i. 
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01 ) 



j+l.i 


+ M • 




a 



where M and T are the moment arm and torque arm found from Equation (10); and 
a p / au , and a T 

j+lfi j+l/i 1 j+1 , i are the elastic axis flexibility vectors 

as found in Equations (9). 

By using the preceding method, systematically for all panels, the complete struc- 
tural influence coefficient matrix is calculated. All panels assigned to the same elastic 
axis endpoint will have the same streamwise rotation. All panels assigned to endpoints 
outboard of the endpoint to which the loaded panel is assigned will have the same rotation 
as the loaded panel . 

The last operation carried out by the structural program is to return the influence 
coefficient matrix to the mainline to be used in the aero-elastic analysis. 

4.2 Method for C omputing the Structural Influence Coefficien t Matrix for Wing-Body- 
Tail Combinations 


For the purpose of computing the structural influence coefficient matrix of a wing- 
body-tail combination, the fuselage is assumed to be represented by a straight torsionally 
rigid beam. The tail surface and wing are represented by an elastic axis as shown in 
Figure 12. The vertical tail is assumed ideally rigid. 

The computer program was developed by adopting the following assumptions: 

t 

(1) The point of intersection of the wing elastic axis and 
the fuselage elastic axis is assumed to be locally rigid. 

See point A of Figure 12. 

(2) Any load on the wing does not affect any port of the fuselage or the 
horizontal tail. This is due to the first assumption. 


(3) The bending deflection of the front fuselage elastic axis 
(i.e. from nose to point A) affects the wing only in a 
region between the wing leading edge and the wing elastic 
axis. The influence of the fuselage angular deflection is assumed 
to decrease linearly from the wing root chord with full deflection 
to a distance at half the structural root chord c s where no influence of 
the fuselage angular deflection is assumed. Thus, as shown in 
Figure 12, only region D of the wing is affected by the front 
fuselage deflection. With the coordinate systems shown in Figure 
12, the angular deflection at any point P' in region D due to 
the front fuselage deflection is given by: 



( 4o ro 0 t> 






( 12 ) 


2 




t. 


Figure 12 Structural Schematic for Wing-Body-Tail Program 


where A ot roQ j. is the angular deflection at the wing root chord 
and is assumed the same as Aa p at a corresponding point P 
(i.e., yp = yp f )on the fuselage elastic axis. In the computer 
program, Aa p is computed by linear interpolation between end- 
points of the fuselage elastic-axis segments. Note that the actual 
deflection at any point on the wing surface is obtained by super- 
imposing the fuselage effect on the wing alone case. 

(4) The same assumption on the influence of the aft fuselage deflection 
applies also to the rear part of the wing surface. See region C of 
Figure 12. 

(5) Any load on the tail surface is first transferred to point B 
and the fuselage deflection is then computed by using beam 
theory. The fuselage deflection affects the horizontal tail 
in two ways. The first is that the angular deflection at 
point B is carried over to the whole tail surface and the 
second is that the tail surface at the junction will follow 
the fuselage in bending. This last effect can be computed 
by assumption (3). For an all-movable tail, this last effect 
does not exist and this is recognized in the program by an 
'if' statement. 

The computer program was written in such a way that efficient and repetitive use 
of the wing-alone structural program can be made. Therefore, the complete structure is 
divided into three parts with symbols , in parentheses: the fuselage (F), wing (W), and the 
horizontal tail (T). If C^py^ represents the angular deflection in the flight direction at 
panel j on the fuselage die to a unit load on panel k of the wing surface, with similar 
definitions for other symbols, then it can be shown that many elements in the structural 
influence coefficient matrix arezero. For example, Cj^py/ = 0 and Cr^yy/ = 0 for all j 
and k. Most elements in Cj^yyp and C:uyyj are also zero. Most elements in C-^pp are 
zero because the front fuselage does not afreet the aft fuselage and vice versa. Tne com- 
plete structural influence coefficient matrix can therefore be expressed as a partitioned 
matrix shown below: 


[C jkFF 1 

[0 ] 

tC jkFT 3 

[C jkWF ] 

[ C jkWW 3 

[C jkWT 

[C jkTF 3 

I 0 ] 

[C jkTT 3 


5. VALIDATION OF THE [C] MATRIX SUBROUTINE 

To test the validity and accuracy of the £cj matrix subroutine, a total of four 
examples were tested. Two of these test cases are presented here. The first, a NACA box 
beam (Reference 2) was used to test the wing-only rC] matrix program. This test case is 
described in Section 5.1 . The second test case, a simple elastic axis model, was used to 
test the wing-body-tail m matrix program and is described in Section 5.2. 



5.1 Wing-Only Test Case 


The thin-walled box beam (Figure 13) presented in Reference 2, was the first 
test wing examined. Figure 14 shows the geometry of the wing and the elastic axis lo- 
cation. Figure 15 shows the El and GJ distributions used. The value of GJ was assumed 
to be infinite in the carry-through section of the box beam, since the mounting of the 
test wing allowed bending but no twisting in this section. A comparison of the computed 
rotations with the experimental results is shown in Figure 16. It is seen that the largest 
discrepancy is of the order of .001 rad (or .0573 degrees). This discrepancy is judged 
to be sufficiently small to consider the program to be working properly. It must be kept 
in mind that it is not possible to state whether the discrepancy is caused by some slight er- 
ror in the computer program or by the beam idealization used in the program. 

5.2 Wing-Body-Tail Case 

No usable test data were available on an actual wing-body-tail configuration. 

For that reason, the computer logic of the extended matrix subroutine was checked 

by performing a computation by hand for a simple wing-body-tail combination as shown 
in Figure 17. The coordinates of control points and elastic axis endpoints are shown in 
Table 1 and 2, respectively. The assumed El and GJ values for each elastic axis segment 
are presented in Table 3. The angular changes in the flight direction (i.e. Aa ) at vari- 
ous locations due to a unit load at some point are tabulated in Table 4. As can be seen 
from Table 4, the computer program produced the same numerical results as computed by 
hand. The computer program is, therefore, considered to be correct. 
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TABLE ; 1 CONTROL POINT COORDINATES FOR FIGURE B. 17 


POINT 

X(in .) 

(m) 

Y(in.) 

(m) 

1 

0.0 

0.0 

0.0 

0.0 

2 

10.0 

0.254 

0.0 

0.0 

3 

20.0 

0.508 

0.0 

0.0 

4 

30.0 

0.762 

0.0 

0.0 

5 

40.0 

1.016 

0.0 

0.0 

6 

50.0 

1.270 

0.0 

0.0 

7 

60.0 

1.524 

0.0 

0.0 

8 

70.0 

1.778 

0.0 

0.0 

9 

80.0 

2.032 

0.0 

0.0 

10 

85.0 

2.159 

0.0 

0.0 

11 

90.0 

2.286 

0.0 

0.0 

12 

26.0 

0.661 

8.5 

0.216 

13 

46.0 

1.168 

9.0 

0.2285 

14 

55.0 

1.396 

18.0 

0.457 

15 

49.66 

1.261 

26.73 

0.679 

16 

54.66 

1.386 

35.39 

0.899 

17 

67.0 

1.701 

6.0 

0.1523 

18 

83.5 

2.120 

10.0 

0.254 

19 

84.303 

2.141 

15.0 

0.381 

Notes: 

Main Wing Structural Root Chord 

- 50.in.(l .270 m.) 




Horizontal Tall Structural Chord = 20. In. (0.508 m.) 
Half-Width of Fuselage at Main Wing = 5. in. (0.127 m.) 
Half-Width of Fuselage at Horizontal Tail = 2. in. (0.0508 m.) 
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TABLE 2 ELASTIC AXIS ENDPOINT COORDINATES FOR FIGURE 


POINT 

X (In.) 

(m) 

Y (in.) 

(m) 

A 

0.0 

0.0 

0.0 

0.0 

B 

10.0 

0.254 

0.0 

0.0 

C 

20.0 

0.508 

0.0 

0.0 

D 

30.0 

0.762 

0.0 

0.0 

E 

40.0 

1.016 

0.0 

0.0 

F 

50.0 

1.270 

0.0 

0.0 

G 

60.0 

1 .524 

0.0 

0.0 

H 

70.0 

1.778 

0.0 

0.0 

I 

80.0 

2.032 

0.0 

0.0 

J 

85.0 

2.159 

0.0 

0.0 

K 

90.0 

2.286 

0.0 

0.0 

L 

40.0 

1.016 

10.0 

0.254 

M 

42.59 

1.082 

19.66. 

0.499 

N 

49.66 

1 .261 

26.73 

0.679 

O 

54.66 

1.388 

35.39 

0.899 

P 

80.0 

2.032 

5.0 

0.127 

Q 

81 .294 

2.065 

9.8295 

0.250 

R 

84.303 

2.141 

13.8225 

0.351 



TABLE 3 ASSUMED El AND GJ VALUES FOR ELASTIC AXIS SEGMENTS FOR FIGURE 17 


SEGMENT 

El (lb- in 2 ) 

Z 

i 

3 

\ G J (lb— i n 2) 

i 

(N-m2) 

A-B 

2.0 x 10 10 

5.74 x 10 7 

l * 

★ 

B-C 

4.0 xlO 10 

1.148 x 10 8 

i 

\ * 

1 

1 

* 

C-D 

5.0 x 10 10 

1.435 x 10 8 

i 

' * 

* 

D-E 

5.0 x 10 10 

1 .435 x 10 8 

i 

! * 

j 

* 

E-F 

5.0 x 10 10 

1.435 x 10 8 

1 

1 * 

| 

* 

F-G 

5.0 x 10 10 

1 .435 x 10 8 

1 

1 * 

i 

* 

G-H 

4.0 x 10 10 

1.148 x 10 8 

i 

: * 

j 

★ 

H-I 

3.0 x 10 10 

8.62 x 10 7 

* 

i 

★ 

I-J 

2.0 x 10 10 

5.74 x 10 7 

t 

* 

* / 

J-K 

1.0 x 10 10 

2.67 x 10 7 

* 

* 

E-L 

1.0 x 10 9 

2.87 x 10 6 

8.0 xlO 8 

2.295 x 10 6 

L-M , 

8.0 x 10 8 

2.295 x JO 6 

6.0 x 10 8 

1 .722 x 10 6 

M-N 

5.0 x 10 8 

1.435 x 10 6 

3.0 x 10 8 

8.62 x 10 5 

N-O 

2.0 x 10 8 

5.74 x 10 5 

CD 

O 

rt 

O 

2.37 x 10 5 

I-P 

3.0 x 10 8 

8.62 x 10 5 

CO 

o 

X 

o 

• 

CM 

5.74 x 10 5 

P-Q 

2.0 x 10 8 

5.74 x 10 5 

1 .0 X 10 8 

2.87 x 10 5 

Q-R 

1.0 x 10 8 

- 2.87 x I0 5 

5.0 x 10 7 

1.435 x 10 4 


Fuselage is assumed torsional ly rigid (GJ = °°) 
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